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Abstract. A multi-component formalism is developed to describe three-body 
systems with nonstatic pairwise interactions and non-nucleonic degrees of freedom. 
The dressed-bag model for NN interaction based on the formation of an intermediate 
six-quark bag dressed by a cr-field is applied to the 3iV system, where it results in a new 
three-body force between the six-quark bag and a third nucleon. Concise variational 
calculations of 3N bound states are carried out in the dressed-bag model including 
the new three-body force. It is shown that this three-body force gives at least half 
the 3N total binding energy, while the weight of non-nucleonic components in the 
and "^He wavefunctions can exceed 10%. The new force model provides a very good 
description of 3A^ bound states with a reasonable magnitude of the aNN coupling 
constant. The model can serve as a natural bridge between dynamical description of 
few-nucleon systems and the very successful Walecka approach to heavy nuclei and 
nuclear matter. 



1. Introduction 

In recent years, few-nucleon theory based on traditional models for the A^A^ force had 
significant success 012]. In particular, a high level of precision for the agreement between 
3A^ and AN calculations and the corresponding experimental data has been reached in 
many cases. Such a degree of agreement obtained without any free parameter could 
promise a description of all properties of few-nucleon systems in the near future by 
solving the exact Faddeev-Yakubovsky equations with the inclusion of modern NN 
potentials and 3A^ forces. 
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However, rather significant disagreements between results of tlie most exact few- 
nucleon calculations and experimental data still remain. Among all such discrepancies, 
we mention here only the most important ones such as the Ay puzzle m. N + d and 
N + ^He scattering [21 E], disagreements at the minima of differential cross sections 
(Sagara puzzle) at £" ~ 150 200 MeV and polarization data ioi N + d N + d, 
N + d^ and N + ^He scattering. The strongest discrepancy between current theories 
and experiments has been found in studies of the short-range NN correlations in the 
^He(e, e'pp) [HI, ^He(7,pp) [Zj, and ^'Re{e,e'NN) [8^ processes. The strong short-range 
correlations found in these experiments are likely due to (at least partially) a three-body 
force acting in a space region where two nucleons are close to each other. Therefore, in 
order to quantitatively describe such correlations, some new types of three-body forces 
should be introduced. 

In addition to these particular problems, there are more fundamental problems in 
the current theory of nuclear forces, e.g., strong discrepancies between the irNN, nNA 
and pNN form factors used in OBE models for the description of elastic and inelastic 
scattering and in the parameterization of 2N and 3A^ forces [HI CHI ^1 Many of 
these difficulties are attributed to a rather poor knowledge of the short-range behaviour 
of nuclear forces. This behaviour was traditionally associated with the vector cu-meson 
exchange. However, the characteristic range of this cj-exchange (for m^^ ~ 780 MeV) is 
equal to about A,^ ~ 0.2 0.3 fm, i.e., is deeply inside the internucleon overlap region. 
Therefore, the quark structure of nucleons should be important at this short range. 
Within six-quark dynamics in its turn, it has long been known fHl CBl ^] that 

the mixing of the completely symmetric s^[6] component with the mixed-symmetry 
s^p^[42] component can determine the structure of the whole short-range interaction 
(in the S-wscve). Assuming a reasonable qq interaction model, many authors (see 
e.g. PHI ^1 1201 1^) have demonstrated that this mixture can result in both strong 
short-range repulsion (associated mainly with the component) and intermediate- 
range attraction (associated mainly with the above mixed-symmetry s^p^ component). 
However, recent studies [201 12I| for A^A^ scattering on the basis of the newly developed 
Goldstone-boson-exchange (GBE) qq interaction have resulted in a purely repulsive NN 
contributions from both s^[6] and s^p^[42] six-quark components. There is no need to 
say that any quark-motivated model for the A^A^ force with vr-exchange between quarks 
inevitably leads to the well-established Yukawa vr-exchange interaction between nucleons 
at long distances. 

Trying to solve the above problems (and to understand more deeply the mechanism 
for the short-range A^A^ interaction), we suggested to add to the conventional Yukawa 
meson-exchange (t-channel) mechanism (at intermediate and short ranges) contributions 
of s-channel graphs describing the formation of a dressed six-quark bag in an 
intermediate state such as \s^ + a) or |s® + 27r)(nil221- It has been shown that, due to the 
change in the symmetry of the six-quark state in the transition from the A^A^ channel to 
the intermediate dressed-bag state, the strong scalar cr-field arises around the symmetric 
6g bag. This intensive cr-field squeezes the bag and increases its density. The high quark 
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density in the symmetric 6g bag enhances meson field fiuctuations around the bag and 
thereby partially restores the chiral symmetry. Therefore, the masses of constituent 
quarks and a mesons decrease. As a result of this phase transition, the dressed bag mass 
decreases considerably (i.e., a large gain in energy arises), which manifests itself as a 
strong effective attraction in the NN channel at intermediate distances. This attraction 
can be described in terms of the OBE model as the effective t-channel cj-exchange. 

However, more accurate calculations of the intermediate-range NN interaction 23] 
within the 27r-exchange model with the vrvr s-wave interaction have revealed that this t- 
channel mechanism cannot give a strong intermediate-range attraction in the A^A^ sector, 
which is necessary for binding of a deuteron and fitting of A^A^ phase shifts. The fact that 
the conventional t-channel 27r-exchange with the vrvr S-wave interaction and reasonably 
soft Attatat cut-off values cannot give the sufficiently strong intermediate-range A^A^ 
attraction has also been corroborated by recent independent calculations p^j. Thus, 
the t-channel mechanism of the a exchange should be replaced by the corresponding 
s-channel mechanism. The contribution of the s-channel mechanism would generally be 
much larger due to resonance-like enhancement | 

Based on this s-channel mechanism, we proposed a new model for the A^A^ 
interaction (referred to as the "dressed bag model" (DBM)) [21123 that provides a quite 
good description of both A^A^ phase shifts up to 1 GeV and the deuteron structure. 
The developed model includes both the conventional t-channel contributions (Yukawa vr 
and 27r-exchanges) at long and intermediate distances and the s-channel contributions 
due to the formation of intermediate dressed-bag states at short distances. The most 
important distinction of such an approach from conventional models for nuclear forces 
is the explicit presence of a non-nucleonic component in the total wavefunction of the 
system, which necessarily implies the presence of new three-body forces (3BF) of several 
kinds in the 3A^ system. These new three-body forces differ from conventionally used 
models for three-body forces. One important aspect of the three-body force should be 
emphasized here. In conventional OBE models, the main contribution to A^A^ attraction 
is due to the t-channel a exchange. However, the three-body force models suggested 
until now (such as Urbana-Illinois or Tucson-Melbourne) are mainly based on the two- 
pion exchange with intermediate A-isobar production, and the cr-exchange either is not 
taken into account or is of little importance in these models. In contrast, cr-exchange 
in our approach dominates in both A^A^ and 3A^ forces; i.e., the general pattern of the 
nuclear interaction appears to be more consistent. 

The aim of this work is just to study both the new type of three-body forces and the 
properties of the 3A^ system with the dressed-bag model (DBM) proposed for A^A^ and 
3A^ forces. In particular, we will show that the role of 3BF changes remarkably compared 
to the conventional OBE models. In the conventional models, the main contribution 
to the 3A^ binding energy comes from the pair A^A^ force, whereas 3BF gives a certain 

X i-channel mechanism can be associated with the direct nuclear reaction where only a few degrees 
of freedom are important, while the s-channel mechanism can be associated with resonance- like (or 
compound nucleus like) nuclear reactions with much larger cross sections at low energies. 
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(although significant) correction of about 15% to the total binding energy. In contrast, 
in our DBM approach the contribution of the new scalar three-body forces turns out to 
be larger by several times and is equal to about half the total 3iV binding energy. 

This paper is organized as follows. In Sect. II, we give a brief description of the DBM 
for the NN system. In Sect. Ill, we develop the multi-component formalism for the 3A^ 
system. In Sect. IV, the general formalism is specified to the one-pole approximation for 
the dressed-bag resolvent, which is used in our calculations. In Sect. V, three types of 
3BF are considered: one-meson exchange (vr and a) between the dressed bag and third 
nucleon and the 2(T-exchange with breaking of the a-loop in the two-nucleon interaction. 
The results of our variational calculation of the 3A^ system are given in Sect. VI and the 
results are discussed in Sect. VII. In the Conclusion we summarize the main results of 
the work. 

2. Dressed bag model for NN forces 

Here, we give a brief description of the two-component dressed-bag model for the A^A^ 
interaction. The detailed description has been presented in our previous papers |^ I22j. 
We consider a system that can exist in two different phase states: the A^A^ phase and 
dressed six-quark bag phase. In the A^A^ channel, the system can be described as two 
nucleons interacting via one-boson exchange; in the Qq + a channel, the system is treated 
as a six-quark bag surrounded by the strong scalar-isoscalar a- field (a "dressed" bag).§ 
Accordingly, the wavefunction of the system consists of two components of different 
nature: 



Here, Hj^p^ includes the direct A^A^ interaction induced by one- and two-pion exchanges 
with the irNN form factors taken with the correct "soft" cut-off parameters 1^ . 
Thus, Hnn describes the peripheral part {at r^N > 1 fm) of the A^A^ interaction. 
H^ff^ is the Hamiltonian of the dressed bag, while the operators Hq^^^mn and H^Nfiqw 
describe the transitions from A^A^ to the {6q+a) channel and vice versa. The Schrodinger 
equation 1 1 reduced to two coupled equations and, by excluding the second component 

§ Full description of the NN interaction at energies E ^ 1 GeV still requires other fields in the bag 
such as 2tt, p and uj but here we employ the version of DBM including only a leading cr-field [22j with 
parameters which are determined either phenomenologically from fit of NN phase shifts (variant (I)) 
or taken directly from the respective cr-loop diagram (the variant (II)). 

II In the force models described here the s-channel mesons are treated relativistically, while the heavy 
Qq bag and nucleons are taken in the nonrelativistic approximation. This is quite reasonable at the 
energy range i?Lab < 0.8 GeV. 




(1) 



The Hamiltonian of the system has the corresponding matrix form 




(2) 
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^Qqj^i one obtains the following equation for the proper NN channel wavefunction with 
the effective Hamiltonian: 

Having obtained the solution of this equation for the "ifNN component, one can also find 
straightforwardly the second component "^sq^: 

^6(^k7 = {E — H^q^) ^ HQg^^NN"^ (4) 

Using the following simple one-pole approximation for the dressed bag resolvent 

where |a) is the 6g part of the wavefunction for the dressed bag and |k) represents the 
plane wave of the a-meson propagation, one obtains an effective interaction in the NN 
channel as a sum of separable terms. ^ We emphasize that the sign of the dressed bag 
wavefunction ^E'eq+o- at low and intermediate energies E < Eq (where Eq is the lowest 
eigenenergy of the dressed bag) is opposite to that of the A^A^ wavefunction. 

To solve equation Q in such an approximation, the knowledge of total Hamiltonian 
iJe-j+cr of the dressed bag, as well as the total transition operator H^N^Qq^, is not 
necessary. Only the projections of the transition operator HNN^Qqj^\a) onto the bag 
states la) are necessary to construct the interaction in the NN channel. These 
projections have been calculated in the microscopic quark-meson model 123- The 
effective interaction V^qN induced in the NN channel by coupling with the intermediate 
dressed bag state is illustrated by the graph in Fig. 1. 




6q 




Figure 1. Effective NN interaction induced by the formation of an intermediate 
dressed bag 



In general, the state index a includes all quantum numbers of the dressed bag, i.e. 
a = {J, M, S, T, Lb, La-}, where Lf,, S, T, J, and M are the orbital angular momentum 
of the 6q bag, its spin, isospin, total angular momentum, and its projection on the z 

% The single-pole approximation for the dressed-bag resolvent should be quite valid at incident energies 
in A''7V-channel Ei^h < 1 GeV because the energy spacing hu> between different excited states of the 
dressed bag (with different parity and number of radial nodes) is expected around 300-350 MeV. So 
that the spacings between excited states in the fixed channel (with the same parity) should be around 
2huj = 600 700 MeV in the cm. system. Hence a few-quantum excitations of the bare 6q bag should 
lie very high in the spectrum. 
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axis, respectively, and is the orbital angular momentum of the a meson. However, 
in the present version of the DBM, only the s-wave state of the 6g bag with the 
configuration is taken into account, so that Lb = 0, J = S , and the isospin of the bag is 
uniquely determined by its spin. The states of the dressed bag with Lg. 7^ should lie 
higher than those with = 0. For this reason, the former states are not included in the 
present version of the model. Therefore, the state index a is specified here by the total 
angular momentum of the bag J and (if necessary) its z projection M: a =^ { J(M)}. 

Thus, the effective interaction in the A^A^ channel V^qN = Hn]^^q^{E — 
Heq^)~^HQ^^NN, after partial- wave decomposition, can be written as a sum of separable 
terms in each partial wave: 

VNqN= Yl ^LL'(r,r'), (6) 

J,L,L' 

with 

Vii'(r, r') = Yl ^f'(r) ^i^'l^) ^l^ir'). (7) 

M 

The energy-dependent coupling constants A^^;(i?) appearing in eq. ((Zj) are calculated 
from the diagram shown in Fig. 1; i.e., they are expressed in terms of the loop integral 
of the product of two transition vertices B and the convolution of two propagators for 
the meson and quark bag with respect to the momentum k: 



JO E-ma-ea{k) 

where e^ik) is the kinetic energy of the dressed bag: ea{k) = k'^/2ma + oj^ik) ~ 
rUa- + fc^/2mo-, rrio- is the reduced mass in the Qq + a channel fh„ = m^rria/ {m„ + rria)'^, 
and uJa{k) being relativistic energy of cj-meson: uJa{k) = ^^/n^^~+~k^. The vertex 
form factors Bj^{k) and the potential form factors ipi^^{v) have been calculated in the 
microscopic quark- meson model I22j. 

When the A^A^-channel wavefunction \l/ArAr is obtained by solving the Schrodinger 
equation with the effective Hamiltonian if at a? -|- V^qN^ the 6gA^-component of the 
wavefunction can be found algebraically from eq. Q: 

\<UE)) = '^i{^.E)\a''^), (9) 

where the function 

^i(^.E).Y.CUE)^^j^^^ (10, 

L 

can be interpreted as the mesonic part of the dressed-bag wavefunction. The coefficients 
Cl{E) in eq. (fTUI) are overlap integrals of the nucleon-nucleon wavefunction and the form 
factors ipjj^^ir) of the separable potential ((Zj): 

Ci{E)= [ ^Mr,E)ifiir)dr. (11) 
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The contribution of the dressed-bag component to the total wavefunction is 
proportional to the norm of '^^q^: 



^6^f = ll«f / ^;;^(k,E)^f^*(k,E)dk 



= Y,Ci{E)Ci;{E)PUE). (12) 

LV 

If the vertex form factors B'l are energy independent, the factors Plj^,{E) are simply 
related to the coupling constants Xj^i^i^E) (see eq.®): 

pHE) = -'^^^. (13) 

The total wavefunction of the bound state \I' must be normalized. Assuming 
that the nucleonic part of the wavefunction "^nn found from the effective Schrodinger 
equation has the standard normalization ||\l/ArAr|| = 1, one obtains that the weight of 
the dressed bag component is equal to 

p = 11^6^11' (14) 



Table 1. Deuteron properties in DBM. 



Model 


^d(MeV) 


Pd{%) 


rm(fm) 


Qd(fm2) 


/id (i^n) 


A5(fm-i/2) 




RSC 


2.22461 


6.47 


1.957 


0.2796 


0.8429 


0.8776 


0.0262 


Moscow 99 


2.22452 


5.52 


1.966 


0.2722 


0.8483 


0.8844 


0.0255 


Bonn 2001 


2.224575 


4.85 


1.966 


0.270 


0.8521 


0.8846 


0.0256 


DBM (var. a) 
Peg = 3.66% 


2.22454 


5.22 


1.9715 


0.2754 


0.8548 


0.8864 


0.0259 


DBM (var. b) 
Peg = 2.5% 


2.22459 


5.31 


1.970 


0.2768 


0.8538 


0.8866 


0.0263 


experiment 


2.224575 




1.971 


0.2859 


0.8574 


0.8846 


0.0263* 



* An average value of the asymptotic mixing parameter r] over a few most accurate results is 
presented here (see refs. [25 1 126 ) l27 1 l2S] ) . 



The model constructed above gives a very good description for the coupled-channel 
^5*1 — ^Di phase shifts, mixing parameter ei and singlet ^5*0 shifts in the energy region 
from zero up to 1 GeV [22]. The deuteron observables obtained in this model without 
any additional parameter are presented in Table 1 in comparison with some other NN 
models and experimental values. The quality of agreement with experimental data for 
the deuteron static properties found with the presented force model, in general, is higher 
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than those for the modem A^A^ potential model such as Bonn, Argonne, etc., especially 
for the asymptotic mixing parameter rj and the deuteron magnetic and quadrupole 
moments. The weight of the dressed bag component in the deuteron is calculated from 
the energy dependence of the coupling constants A^^, (see eqs. (fT^ and (fT^ and is 
varied from 2.5% to 3.6% in different versions of the model [01 122] ■ 



3. DBM for the three-nucleon system 



When generalizing the above model to the 3A^ system, one meets two difficulties: the 
appearance of the energy dependent coupling constants \j^]^,{E) and the presence of 
non-nucleonic degrees of freedom in an explicit form. Thus, we suggest that the 3A^ 
system can be found in the five different states (in three phases): NNN, Ni + {6q^^^^a), 
N2 + (6g(3^V), A^s + (6g(^^V), and 9q + a. Therefore, the total wave function can be 
represented as the five-row Fock column 



(2) 
6gN 
(3) 



(15) 



The nucleonic component \l/37v here describes the system as three nucleons 
interacting at large distances via (one- and two-) boson exchanges. In the (GgA^)'-*-' 
channels, there are the dressed Qq bag and ith nucleon [i = 1,2,3). Hereinafter, the 
symbol a in the notation of these channels is omitted for brevity. The wavefunction 
component \I/gg^ depends on the dressed bag variables, including the momentum of the 
(j-meson, and the 6q—N relative motion variable (the Jacobi coordinate or momentum). 
The fifth component of the total wavefunction ()15|1 describes a dressed nine-quark bag. 

The admixture of the dressed 9g-bag in 3A^ system and in "normal" nuclei should 
be rather low, at least one order of magnitude lower as compared to the weight of the 
6q bag. This can be easily explained by the fact the 9g-clusters in nuclei can be only 
generated from 6q and 3g-clusters at their fusion which has rather low probability due 
to very high kinetic energy of 9q system. So the weight of 9q configuration must be a 
small fraction of that for Qq + 3q component. Thus one can expect only a very minor 
effects of this 9q configurations except the region of very high momentum transfer. 

Therefore, we neglect this component in present study and consider below the 4- 
component system with the (4x4)-matrix Hamiltonian 



H 











\ 




Hi 










H2,3N 





H2 







\ H3 3N 








H3 


/ 



(16) 



Here, we use the brief notation for the following parts of the Hamiltonian: the operator 



H- 



3N,i 



= H. 



describes the transition from the {{6q^^''^ + a) + Ni} channel to 
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the pure three-nucleon channel; the Hamiltonian Hi = H^q^^- describes the interacting 
system of the dressed 6q bag and ith nucleon. 

In our 3A^ dynamics, we neglect the direct transitions between different 6q — N 
channels. Such channels are coupled only via an intermediate 3N configurations. 
This assumption looks absolutely natural, while the incorporation of direct transitions 
between such channels requires some exotic new mechanisms for the coupling. 

The Hamiltonian H^j^i of the nucleonic channel includes the kinetic energy of the 
relative motion of three nucleons and the direct A^A^ interactions (OPE and TPE with 
soft cut-offs) that enter into the two-body Hamiltonian: 

3 

hNN=t + VMN H3N = T + J2yNN- (17) 

i=l 

Hereinafter, we use lowercase letters {g, h, etc.) for two-body operators and capital 
letters for three-body ones. 

The (6gA^ 3A^)-transition operators H^n,! do not affect the third nucleon and 
therefore take the form: 

H3N,i = hfl,^ ® 1«, where {tjk) = (123), (231), (312), (18) 

and h%^^ 

6,jfCT the two-body DBM transition operators for pair {jk). 
Since direct transitions between different QqN channels are absent, the exclusion of 
the QqN components is completely similar to the two-body case: hence one can rewrite 
the 4-component Hamiltonian in the two-component form 

H={ l'"" I'''''''' ] , (19) 

where H^gN is the diagonal (3x3)-matrix operator with the elements Hi. The resolvent 
of i^egAf is also a diagonal (3x3)-matrix. 

The Hamiltonian Hi describes the dressed bag Qq{jk) (formed by the fusion of the 
jth and kth nucleons) interacting as a whole with the third (zth) nucleon via exchanges 
by mesons (vr, a). Just this interaction results in the new three-body force arising in 
the two-phase model. The Hamiltonian 

Hi = he^{jk) © /i6g,7v(«) , (20) 

is the direct sum of the dressed-bag Hamiltonian hQq^{jk) and single-particle 
Hamiltonian of the ith nucleon h^gj^a), which includes its kinetic energy and the 
potential energy of its interaction with the bag: 

hq,N(') =U + VQg^p^(i) (21) 

It depends only on the single Jacobi coordinate (or momentum) of the third (ith) 
nucleon. Therefore, the resolvent of Hi is equal to the convolution of two subresolvents: 
the first for the dressed-bag Hamiltonian {E — hQg^{jk))~^ and the second for the single- 
particle Hamiltonian ()2H1 : 



1 f°° 

{E-Hir' = -— / (e~heg^{jk)r'g,g,MAE-e)de, 



(22) 
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where 

g,,,^XE) = {E~h,^^^,.,)-\ (23) 

Excluding the QqN channels from the Schrodinger equation for the total four- 
component wave function, one arrives at the following effective Hamiltonian T^sa? acting 
in the 3iV channel: 

3 

'H'iN = + ^ H^N,i{E — Hi)^^Hi^3N (24) 
1=1 

which has the form of the standard three-body Hamiltonian with pairwise energy- 
dependent interactions: 

3 

Ti-SN = T + ^(V^*]^ + y3N,6qN,3N)'^ (25) 
i=l 

This form of the effective Hamiltonian is quite suitable for the standard Faddeev 
reduction. Although the effective interaction is the sum of three "pairwise" interactions, 
each term V^^^^^j^^j^ describes not a genuine two-body force but a three-body force 
including the interaction between the Qq bag and third nucleon. Moreover, even if this 
interaction is disregarded, the potential Kj^^g^^yg^ still depends on the momentum of 
the third nucleon. This dependence on the momentum of the third particle reduces the 
strength of the effective interaction between two other particles due to a specific energy 
dependence of the coupling constants. Actually, the form of the effective Hamiltonian 
()24|) implies that there are no pure two-body forces in the 3A^ system in our approach, 
except the peripheral part of the OBE interaction. 

After solving the Faddeev equations (or the Schrodinger equation for bound states) 
with the effective Hamiltonian ()24|) for the 3A^ wavefunction, the "bag" components of 
the total wave function are "recovered" by means of the relationship: 

= iE- H,)-'H,,sN^sN. (26) 

4. Single-pole approximation for the dressed bag resolvent 

The momentum representation is more appropriate for description of the 3A^ system in 
the case of DBM. We will employ the same notation for functions both in the coordinate 
and momentum representations. In the single-pole approximation for the dressed-bag 
resolvent, the resolvent {E — Hi)~^ fl22|) can be reduced to a sum of factorized terms: 



{E-H,)-' = Y, j |a('U)6;6,,^^(E-m„-£.(A;))(a«,k|dk 



(27) 



and the effective 3A^ interaction V^^^ 3^ becomes the sum of integral operators with 
factorized kernels: 

^3S,6.Ar,3iv(Pi, p'i, qi> q'i; e)= Y1 ^i^'ipOwUliqu 4; i?)^i'^'(p:), (28) 

JM,J'M',L,L' 
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where Wlf, is expressed in terms of the integral of the product of the vertex functions 
and one-particle resolvent (PHj) : 

(29) 



If the interaction between the 6q bag and third nucleon is neglected, the resolvent Qeq^N 
becomes the simple propagator: 

1 m]\fma 



m 



-. (30) 



E — ma — ea{k)—q^/2m^ rriN + nia 

Then, the effective interaction reduces to the sum of two-body separable potentials with 
the coupling constants depending on the total three-body energy E and third-particle 
momentum qf. 

2m' 



v^i5,6,;v,3iv(Pi,pUi,q'i;^) = 5(qi-q:) E ^f'(Pi) ^iL'(^ - 7^)^r(p'i), (3i) 



J,M,L,L' 

2 

where the coupling constants X^^, are determined by eq.® with replacing E hj E — 
In other words, the effective interaction in the three-body system reduces to a sum of 
three effective pairwise interactions V^gN depending on the total three-body energy and 
momentum of the third (spectator) particle: 

3 2 

VsN,6,N,3N E5(q,-q',)l^i;)^(pi,p'i;E-|^). (32) 



i=l 



It should be noted that this result does not depend on specific choice of interaction 
model, but is quite general provided that there are not direct transition between different 
6qN channels and the transition operators take the form (fTHj) 



a 



1 1 OPE + TPE 




""f •k; 






, 6q ^"^^ 








• 










|3BF 






* 



1 



Figure 2. Different interactions in the 3A^ system for one of tliree possible 
combinations (1+23): tlie peripheral two-nucleon interaction is due to OPE + TPE, 
effective two-body interaction vj^^jj^ is induced by the formation of dressed six-quark 
bag and meson-exchange three-body force (3BF). 



When using the effective interaction (jH^ . one must also include an additional three- 
body force due to the meson-exchange interaction between the dressed bag and third 
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nucleoli. The pattern of different interactions arising in the 3A^ system in such a way is 
illustrated in Fig. 2. 

In the single-pole approximation, the components of the "bag" wave function ()2fij) 
can be expressed in terms of the nucleonic component as 

*g.(k.q,;£)= E |a-'.''.>-?iM^p2iL^, (33) 

where X'l'^'i^i) ^ire the overlap integrals of the three-nucleon component of the 
wavefunction and the potential form factors v^^'.^*: 



xir' (q^) = J (p^) ^3^(Pi' qi) dpi- (34) 

These overlap functions depend on the momentum (or coordinate), spin, and isospin 
of the third nucleon. For brevity, the spin-isospin parts of the overlap functions and 
corresponding quantum numbers are omitted unless they are needed. In eqs.(|HH |) -(|H ^ 
and below, we keep the index i in the quantum numbers L and J in order to distinguish 
the orbital and total angular momenta in the form factors from the respective angular 
momenta of the whole 3A^ system. 

It should be noted that the angular part of the function X^l\s\i) i'^ eq. p^ is not 
equal to 1^-^,/^ (^)- This part includes other angular momenta due to coupling between 
the angular momenta and spins of the dressed bag and third nucleon. The total overlap 
function (without isospin part) can be written, for example, as (we use here the same 
letter for notation of both the total wavefunction and its radial part): 

xif''(qO= E xt,(g.)U"^^^.M,|jM)3;,^7/(q,). (35) 

Here, J and M are the total angular momentum of the 3A^ system and its z projection, 
respectively; and J are the orbital and total angular momenta of the third (ith) 
nucleon, respectively; Jj is the total angular momentum of the bag (Jj = S'j, because in 
our case the bag has zero orbital momentum); and Li is the orbital angular momentum 
of a nucleon pair (j/c) related to the transition vertex [NN — >■ Gg-t-a). For the total 
angular momentum J = 1/2 (the ground state of and ^He), the orbital angular 
momentum of the third nucleon (in the c.m.s. of the entire system) can take the 
values Jj ± 1 (we mean here the 6giV component only). Therefore, li in our present 
model can be equal to or 2. 

The total norm of all three 6gA^ components for the 3A^ bound state is determined 
by the integral 

\\^e,Nr = ^J2\\^''rJ2 / / -— dk UA(q.)dq.(36) 

If the vertex functions Bj^ik) are energy independent (it was just conjectured when 
deriving the effective Hamiltonian P^ ). the internal loop integral with respect to k in 
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eq.(jnSl) (in braces) can be replaced by the energy derivative of A;^^,: 
5£(k)i?^Uk) ^, _ d /• _Bi^(kK^(k) 



(E - - ^ - e„{k)f J E~m^-^- e^{k) 

= -4-Xi\,iE--^). (37) 
dE ^'^^^ 2m' ^ ' 

Thus, the weight of the QqN component in the 3A^ system is determined by the same 

energy dependence of the couphng constants Xiui^) the contribution of the 6q 

component in the NN system but at a shifted energy. 

Using eq. ()37|) . the norm of 6qN component can be rewritten eventually as 

ll*65iv|r = 3^||a^'||^ J xi\{q^(^-^\iiE-q^/2m)^xi\{cLi)dci,, (38) 

Due to the explicit presence of the meson variables in our approach, it is generally 
impossible to define the wavefunction of the third nucleon in the 6qN channel. However, 
by integrating \E'6gAf(k, q) with respect to the meson momentum k, one can obtain the 
momentum distribution of the third nucleon in the 6qN channel weighted with the a- 
meson momentum distribution. Based on eq. (j38|l . we can attribute the meaning of the 
wavefunction of the third nucleon in the 6qN channel to the quantity 



^.''•(q.) = y (-^Ai;,,(E-g,V2m))xif''(q.), (39) 

With this "wavefunction" , one can calculate the mean value of any operator depending 
on the momentum (or coordinate) of the third nucleon. We note that the derivative 
—dX/dE is always positive. 



5. Three-body forces in the dressed bag model 

In this work, we employ the effective interaction ()3H1 and take into account the 
interaction between the dressed bag and third nucleon as an additional three-body force 
(3BF). We consider here three types of 3BF: one- meson exchange (vr and cr) between the 
dressed bag and third nucleon (see Figs. 3a and 3b) and the exchange by two a-mesons, 
where the third-nucleon propagator breaks the cr-loop of the two-body force (Fig.3c). 




Figure 3. The graphs corresponding to three new types of three-body force 
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All these forces can be represented as some integral operators with factorized kernels 
similar to the effective interaction V^^^ggArsAf 

^^^^(^^(Pi, P- qi, q^; E)= ^i^(pO ^^^iy//;(qi, q^; E) ^f,^^'!?;), (40) 

JM,J'M',L,L' 

Therefore, matrix elements for 3BF include only the overlap functions, and thus 
the contribution of 3BF is proportional to the weight of the QqN component in the total 
wavefunction. To our knowledge, the first calculation of the 3BF contribution induced 
by one-pion exchange between the 6g bag and third nucleon was done by Fasano and 
Lee |2ni in the hybrid quark-compound bag model (QCB) using perturbation theory. 
They used the model where the weight of the 6g component in a deuteron is ca. 1.7%, 
and thus they obtained a very small value of -0.041 MeV for the 3BF OPE contribution 
to the 3A^ binding energy. Our results for the OPE 3BF agree with the results obtained 
by Fasano and Lee, because the OPE contribution to 3BF is proportional to the weight 
of the 6g component, and it should be in our case at least twice as large as in their 
calculation. However, we found that a much larger contribution comes from scalar- 
meson exchanges (OSE and TSE). We emphasize that, due to (proposed) restoration 
of chiral symmetry in our approach, the cr-meson mass is ca. 400 MeV, and thus the 
effective radius of the cr-exchange interaction is not so small as that in conventional OBE 
models. Therefore, we cannot use the perturbation theory anymore to estimate the 3BF 
contribution and have to do the full calculation including 3BF in the total three-body 
Hamiltonian. 



5.1. One-meson exchange between the dressed bag and third nucleon 

For the one-meson exchange (OME) term, the three-body interaction '^^^Wlf^, takes 
the form: 

^^^iy//;(q,q:;E) = 

dk ^ W(qi, qO ^ . (41) 

E - nia- 1^- ea{k) E - nia - ^ - e„{k) 

Therefore, the matrix element for one-meson exchange can be expressed in terms of the 
wave function of the "bag" components ^E^ggV- 

(vI/3^|OME|vI/3^) = m,,^\V'''^''\%,^) (42) 

Including spin-isospin variables (it is important for evaluation of pion-exchange 
interaction), we can write down ^'^^^^ as 

^» - V '(q) \T, I -TT) (A3) 

JiM.M ^ ~ ^''^'^) 

where is the bag isospin, T and are the total isospin of the system (T = 1/2 for 
and ^He) and its 2;-projection respectively and Xl,^^'{^) the overlap functions (j^Sl)- 
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When calculating the matrix elements for OME, one can use a similar trick as 
in the calculation of the norm for the QqN component. It enables us to exclude the 
vertex functions Bj^^ (k) from the formulas for matrix elements. Replacing the product 
of propagators in the integral with respect to the meson momentum k (eq. ()41|) ) by their 
difference, one obtains the following expression free of the vertex functions: 



/ 



= AAl^,(g,g') (44) 

This quantity is the finite-difference analogue of the derivative of A with respect to q^. 
In the present calculations, we employed a rational approximation for the energy 



dependence of A/ 



Ai..(i^) = Ai,,(0)§±^, (45) 

where the parameters Eq and a are taken to be the same for all A's. We found that this 
simple rational form can reproduce quite accurately the exact energy dependence of the 
coupling constants \'[j^i{E) calculated from the loop diagram in Fig. 1. For an energy 
dependence such as in eq.(H3|). the energy derivative and AA(g, g') have the form 

-AA(B) = A(0)Ml+a)^^^; (46) 

AA(,.,0 = A(O)Bo(l + a)^3^Jf37v5;i- ^ 

Thus, the matrix elements for OME can be found without using the vertex functions 
5/(k) explicitly. 

5.1.1. One-pion exchange 

For one-pion exchange, we take the interaction operator in the standard form 

= -(iS]l(-"'P)^T^(S^P)(-'"T.). P = 1'. (48) 

where o"*^*-* and r*^*^ are the spin and isospin variables of the third (ith) nucleon, whereas 
Sd and Td are the operators of the total spin and isospin of the 6q bag, respectively. 
We found that the contribution of OPE is so small that it is sufficient to include here 
only S waves. In this case, only the central part of the OPE interaction remains: 

= ^7^^^7^^(cr«Sd)(r«Td)^^. (49) 

The spin-isospin matrix element is nonzero only for a singlet-triplet transition: 

(5,=o,r,=i|i(c7«Sd)(T«Td)|Sd=i,rd=o) = ^, (50) 
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Then, the matrix element of the OPE contribution for s waves takes the form 



(OPE), = -fl^^.J\l,mU^) Eo{l + a) 

Here, we take the vertex functions Bq and Bq differing from each other only by a 
constant. Therefore, using eg . (j44j) . one can exclude these functions from the formula 
for the matrix elements. 

We did not introduce here any cut-off factor for OPE, because the overlap functions 
truncate the OPE interaction at large q values. It should be noted that the overlap 
functions Xo ^^^d Xo singlet and triplet bags, respectively, are very similar in shape 
and magnitude but have opposite signs (see Fig. 5). Therefore, the OPE contribution 
(jnH) is always negative. In our case, it is equal to only —0.1 MeV and thus generally 
agrees with the result obtained in ref . ||23 when the latter is rescaled to the larger weight 
of the 6g component in our case. 

5.1.2. One-sigma exchange (OSE) 

The scalar meson exchange operator does not include any spin-isospin variables. 
Therefore, eq. ()41|) for W///' can be simplified and, in view of the energy dependence 
given in eq. (jlSj) . reduces to the form 

^^^;y//;(qi, q^; E) = 5jj.Ai^,(0) Eo{l + a) 

V ^ -91nN ^ 

Thus, the matrix element of OSE 3BF is determined by the formula 
{0SE) = -3gl^^ XL'WEoil + a) 

JM,L,U 

In the actual calculations, we used the "light" cr-meson mass m„ = 400 MeV. 

The physical meaning of the OSE three-body force is easily understood: it 
represents the interaction of the third nucleon with the a-meson cloud of the bag. In 
view of the enhancement of the cx-field around the dense 6q bag, this contribution is 
very important for our understanding of the properties of 3A^ and heavier nuclei. 

5.2. Two-sigma exchange (TSE) 

The two-sigma process shown in Fig.4 also contributes significantly to 3BF. This 3A^ 
interaction seems less important than the OSE force, because this interaction imposes 
a specific kinematic restriction on the 3A^ configuration."*" 

"•" It follows from the intuitive picture of this interaction that this force can be large only if the 
momentum of the third nucleon is almost opposite to the momentum of the emitted cr-meson. Thus, 
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Figure 4. Three-body force due to two-sigma exchange 



The operator of this interaction includes the vertex functions for the (A^A^ 6g-|-cr) 
and reverse transitions so that these vertices cannot be excluded similarly to the case 
of OME. Therefore, we have to choose some form of these functions. It is naturally 
to suppose that these vertices are the same as those assumed in two-body DBM; i.e., 
they can be normalized by means of the coupling constants A(-E), which, in turn, are 
chosen in the two-nucleon sector to accurately describe A^A^ phase shifts and deuteron 
properties (see eq. ()56|) for vertex normalization). We use the Gaussian form factor for 
these vertices: 

5/(k) = Eq^^ ^ (54) 

where k is the meson momentum and the parameter b is taken from the microscopical 
quark model [22]: 

= ^bl; bo = 0.5 fm. (55) 
Then, the vertex constants Bq should be found from the equation: 

dk ,^ ° ° ' ^,. =XiAE), (56) 



where X'li^^E) are the coupling constants employed in the construction of the DBM in 
the A^A^ sector and are fixed by A^A^ phase shifts. 

For the aNN vertices, we also take the Gaussian form factor: 

QaNNe-'"'''', with = hi (57) 



Then, the box diagram in Fig. 4 can be expressed in terms of the integral with respect 
to the momentum qo of the third nucleon in the intermediate state: 

TSEtjtJJ' I. rp\ _ ^ „2 TjJL T^JL' 

^LU 1^! q , -C^J - J J' 9 aNN ^0 ^0 

(2^0^ J ^"^V^ + (q„_q)2 E _A ml + (qo-q')^' ^^^^ 

" 2 m 

a specific 3A^ kinematic configuration is required when two nucleons approach close to each other to 
form a bag, while the third nucleon has a specific space localization and momentum. 
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Thus, the matrix element for the total contribution of TSE takes the form 

(TSE) = 3 /xf'(q)^^^W^//'(q,q';i^^)xii''(q')dqdq'. (59) 

JM,L,L' 

After the partial wave decomposition, these six-dimensional integrals can be reduce to 
two-dimensional integrals, which are computed numerically by means of the Gaussian 
quadratures. 

6. Results 

The present calculations of 3A^ bound states have been made with the two variants 
of the DBM, versions I and II, one of which being used in our previous two-nucleon 
studies [22] • For the other variant, we use a weaker energy dependence X{E), which is 
taken directly from the loop integral (eq. ()5fij) ) corresponding to the diagram in Fig. 1. 
The effective three-nucleon Hamiltonian in the 3A^ channel includes both two-body and 
three-body forces (although, as was noted above, such a division is justified only for the 
Yukawa (t-channel) exchange contribution): 

i=l i=l 

The parameters of two-body potentials V^^^ (for the two versions) provide a quite 
good description of the deuteron observables A^A^ and phase shifts in the wide energy 
interval 0-1 GeV. The parameters of the three-body force V^^^^ were discussed in 
the previous section. We took into account both the OSE and TSE contributions to 
3BF. The contribution of OPE interaction between the dressed bag and third nucleon 
to the 3A^ binding energy was found to be only about -0.1 MeV. Therefore, we did 
not include this term in the complete variational calculation but estimated it with 
the perturbation theory (PT). The present 3A^ calculations were made using a highly 
accurate variational method and the antisymmetrized Gaussian basis [HHI l!?T] with non- 
linear scale parameters. This basis is extremely flexible and has an enormous advantage 
for such tedious calculations. It enables us to obtain completely analytical formulas for 
the matrix elements of effective two- and three-body forces, except two-sigma exchange 
(these formulas will be presented in another paper). In the present calculations, we used 
the variational basis including even partial waves with the total orbital angular momenta 
L = 0, 2 and took into account the A^A^ force components for the ^5*0 and ^5*1 — ^Di 
channels. The dimension of the basis was increased until the converged results were 
reached. The results presented in Table 2 are obtained with the basis that includes six 
three-body partial components and has the total dimension A^ = 688. 

Table 2 presents our results for the following static properties of ^H: binding energy 
E, D-wave percentage Pd, and weight of the dressed "bag" 6qN component P^gN- In 
addition, this table shows the contributions of the following individual parts of the 
Hamiltonian to the total three-body expectation value: three-nucleon kinetic energy 
T, two-body effective force V^I^, and three-body force due to one-sigma (Vjf^) and 
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Table 2. Results of 3A^ calculations with two- and three-body forces for two versions 
(I and II) of the dressed-bag model. 





E 
MeV 


Pd % 




Individual contributions to H, MeV 


T 


-i/NN 
^ NqN 


^OSE 


T/3BF 
^TSE 




Version I, g^NN = 9.5 


3H(i) 


-5.83 


6.44 


7.87 


75.0 


-64.45 


-1.54 (PT)* 


-0.43 (PT)* 




3h(2) 


-4.14 


5.81 


6.24 


58.1 


-49.88 


-0.83 (PT)* 


-0.22 (PT)* 


-0.01 


3h(3) 


-8.326 


6.84 


10.84 


110.7 


-92.9 


-4.72 


-2.03 


-0.1 


3He(3) 


-7.588 


6.80 


10.66 


108.2 


-91.0 


-4.55 


-1.99 


-0.1 


Version II, Qunn = 8.6 


3H(i) 


-6.12 


6.67 


5.45 


79.0 


-67.82 


-1.30 (PT)* 


-0.56 (PT)* 




3h(2) 


-5.01 


6.23 


4.67 


66.9 


-57.3 


-0.87 (PT)* 


-0.36 (PT)* 


-0.01 


3^(3) 


-8.358 


7.06 


7.31 


110.8 


-94.2 


-2.79 


-1.66 


-0.1 


3He(3) 


-7.565 


7.00 


7.14 


107.8 


-91.8 


-2.66 


-1.58 


-0.1 



- without 3BF, without g^-dependence, 

(2) - without 3BF but with q^-de\ieiidi.ei\ce^ 

(3) - with 3BF and with g^-dependence, 



* These values have been found with the 3A^ wavefunctions in those calculation the 3BF 
have been omitted, i.e. within PT. The comparison of these entries with those found in the 
complete calculation (see rows denoted ^H^'^^) shows with evidence inapplicability of the PT. 

two-sigma exchanges iV^g^). The contribution of 3BF due to OPE calculated with 
perturbation theory (PT) is also given in Table 2. The other details of the results 
presented in the Table, especially related to the ^He properties, are discussed in the 
forthcoming paper. 

For both variants of the model, we present also the result calculated disregarding 
both 3BF and the dependence of the effective two-body force V^qN on the momentum 
of the third nucleon (see the first and fifth rows of Table 2). The results in the second and 
sixth rows of Table 2 are obtained including the dependence of V^^gA^, but disregarding 
3BF. The percentages of the D-wave Pd and QqN component P%qN given in Table 2 were 
obtained with incorporation of the three 6giV components; i.e., these values correspond 
to the normalization of the total (four-component) wavefunction of the system to unity. 

7. Discussion 

The 3iV results presented in the previous section differ significantly both from the 3A^ 
results found with the conventional models for A^A^ and 3A^ forces (based on Yukawa's 
meson exchange mechanism) and from the results obtained in the framework of hybrid 
models , which include the two-component representation of the NN wavefunction 
^ = "^NN + ^6(j- It is convenient to discuss these differences in the following order. 

(i) We found that the dependence of pair A^A^ forces on the momentum of the third 
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particle in the 3A^ system is more pronounced in our case than in other hybrid 
models [SH EHl OH ES]: the 3A^ binding energy decreases by more than 1.7 MeV, 
from 5.85 to 4.14 MeV (cf. the first and second rows in Table 2). From the more 
general point of view, it means that, in our approach, pair NN interactions (except 
Yukawa OPE and TPE terms), being "embedded" into a many-body system, loose 
their two-particle character and become substantially many-body forces (depending 
on the momenta of other particles of the system). 

(ii) Due to such a strong dependence (of "repulsive" character), the 3A^ system 
calculated including only pair forces turns out to be strongly underbound {E = 
—4.14 MeV). In other words, the"pairwise" A^A^ forces (including their 
dependence on the momenta of the third nucleon) give only about half the total 
3A^ binding energy, leaving the second half for the 3BF contribution. Therefore, 
the following question is decisively important: can the three-body force (inevitably 
arising in our approach) give the large missing contribution to the 3A^ binding 
energy? Usefulness of the developed model for the description of nuclear systems 
depends directly on the answer to this important question. It is appropriate here 
to remind that in the conventional 3BF models such as Urbana-Illinois or Tucson- 
Melbourne, the contribution of 3BF to the total 3A^ binding energy does not exceed 
1 MeV; i.e., this contribution can be considered as some correction (~ 15%), 
although it is significant for the precise description of the 3A^ system. 

(iii) Fortunately, the contribution of 3BF induced by cr and 2a exchanges enables one 
to fill this 4.4 MeV gap between the two-body force contribution and experimental 
value. In fact, including both one- and two-sigma exchange contributions to 3BF, 
taken with the same parameters as in the initial NN-force model and a quite 
reasonable coupling constant gaNN = 8 -h 9.5, we obtain the 3A^ binding energy 
that is very close to the experimental value (see rows 3,4 and 7,8 in Table 2). Thus, 
the presented force model leads to a very reasonable binding energy for the three- 
nucleon system but with the strongly enhanced (as compared to the traditional 3A^ 
force model) contribution of three-body forces. 

(iv) The contributions of pair and different three-body forces to the total 3A^ binding 
energy for are given in the third and fourth rows of Table 2. From the results 
presented in these rows, one can conclude that the total 3BF contribution to the 3A^ 
binding energy dominates and, in fact, determines the structure of the '^H ground 
state*. Moreover, comparing the third and fourth rows of the table, one can see a 
"nonlinear" effect of self-strengthening for the 3BF contribution. The comparison of 
the results presented in the third and fourth rows of the table (see the second, fourth, 
and fifth columns) shows clearly that the binding energy is almost proportional to 
the weight P^qN of the 6qN component in the total 3A^ wavefunction. Thus, when 
the weight of the 6qN component increases, the 3BF contribution, which is related 

* It should be noted that the relative contribution of the pair effective force YnqN to the binding energy 
decreases noticeably when including 3BF (due to strengthening the dependence of pair forces). 
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directly only to this component of the total wavefunction, increases accordingly. 
However, the enhancement of the pure attractive 3BF contribution squeezes the 3A^ 
system and thus reduces its rms radius, i.e., the mean distance between nucleons, 
which, in turn, again increases the weight of the QqN component. In other words, a 
some chain process strengthening the attraction in the system arises. This process 
is balanced both by the weakening of the effective pair interaction due to the 
dependence and by the repulsive effect of the orthogonalizing pseudopotentials 
included in each pair interaction. 

There are two another important stabilizing factors weakening the strong three- 
body attraction in the 3A^ system. First, the generation of the short-range repulsive 
vector cij-field, where all three nucleons are close to each other ^j. Since the uj- 
meson is heavy, this field is located in the deep overlap region of all three nucleons. 
In this study, we omitted the three-body contribution of this repulsive ti;-field. This 
repulsive contribution will keep the whole system from the further collapse due to 
the strong attractive 3iV force induced by the scalar field. 

The second factor slightly weakening the effective 3A^ attraction is associated 
with the conservation of the number of scalar mesons generated in the 2N and 3A^ 
interaction process. The problem is that 2(T-exchange giving the 3BF contribution 
(see Fig. 4) arises due to the break of the cr-meson loop, which induces the main 
2N force. In other words, the a-meson generated in the transition of pair nucleons 
from the NN phase state to the 6g state is absorbed either in the 6g bag with 
closing the loop or by the third nucleon resulting in the 3BF contribution. Thus, 
the appearance of such a 3BF should weaken attraction between nucleons in the 
pair. We carefully estimated the effect of the meson-number conservation for the 
TSE contribution on the total 3A^ binding energy. Its magnitude occurred to 
be rather moderate in the absolute energy scale (ca. 0.3 - 0.4 MeV) but quite 
noticeable within the whole TSE contribution. However, when the total nucleon 
density increases (and the relative TSE contribution also increases), the effect is 
enhanced. 

(v) Dependence of the two-body coupling constants \{e) upon the average momentum 
of other nucleon in three-nucleon system (see e.g. eq. 1)311) can be interpreted 
generally as a density dependence of the two-body force in many-nucleon system. 
Thus, one fixes (e.g. at energy e = Ed) the value of energy dependent coupling 
constant \{E—q^ /2m) of our two-body force, i.e., if one disregards its g^-dependence 
(this g^-dependence leads to a weakening of the two body force in a many-nuclear 
system when q^ is rising), than the neglected g^-dependence must be compensated 
by an additional repulsive density- dependent effective three-body force. 

On the other hand, it is well known from the Skirme-model calculations 
of nuclei that just similar repulsive phenomenological density-dependent three- 
body force should be added to conventional 2N- and 3A^-forces to guarantee the 
saturation properties of heavy nuclei. Thus, in this respect the present force 
model is also in a qualitative agreement with phenomenological picture of nuclear 
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0.18 




-0.04 



Figure 5. QqN component of the wavefunction: the triplet and singlet s-wave 
overlap functions Xq{p) (solid line) and Xo(p) (dashed line) and the corresponding 
"quasi- wavefunctions" ipoip) (dash-dotted line) and V'o(p) (dot-dot-dashed line). The 
triplet functions are multiplied by a factor of (-1). 



interactions. 

(vi) Figure 5 shows two types of the s-wave wavefunctions in the 6qN channeh 
the triplet xlip) (sohd hne) and singlet Xo(p) (dashed line) overlap functions 
given by eqs. (jHH ESI) ^^e corresponding "quasi- wavefunctions" ipQ^p) and ipQ^p) 
normalized according to eq. (jSHI)- The triplet functions are multiplied by (-1) for 
convenience. It is evident that there is an inner node (at r„ ~ 0.56 fm) in all 
s-wave components of the quasi-wavefunctions and the overlap functions. These 
inner nodes are stationary when the total energy increases and the Nd scattering 
problem is considered. As is seen in Fig. 5, the overlap function Xoip) fo^' the triplet 
QqN channel has a more extended tail than the singlet component Xoip)- This is 
due to the fact that the 3A^ system in the spin-triplet channel (coinciding with 
deuteron channel) is less bound, by about 2 MeV, than in the spin-singlet channel. 
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8. Conclusion 

In this paper, we have developed a formalism for a multi-component description of 
the three-nucleon system within the new approach to the 2N and 3iV interactions 
based on the dressed dibaryon and cr-field generation. It has been shown that the 
DBM applied to the 3A^ system results automatically in a new three-body scalar force 
due to the interaction between the dressed dibaryon and third nucleon. This force 
plays a crucial role in the structure of few-nucleon systems. Our accurate variational 
calculations have demonstrated that new 3BF gives half the 3N binding energy, whereas 
the 3BF contribution in the traditional NN-iorce approaches gives about 15% of the 
total binding energy. Thus, the suggested approach to the NN interaction can lead to 
significant revision of relative contributions of two- and many-body forces in nuclear 
systems. 

It is crucially important that the DBM gives an 8 - 11% non-nucleonic component 
in the 3N wavefunction, while this component in the deuteron is equal to only about 
3%, leading to a reformulation of many effects in few-nucleon systems and other nuclei 
as well. It is probable that the weight of such non-nucleonic components in heavy nuclei 
can be even higher with an increase in the mass number and nuclear density. 

Generalization to 4N (and more) systems can be done rather similarly to the 
case of conventional force model. Having 2N- and 3A^-forces fixed from 2N- and 3A^- 
calculations one can straightforwardly consider systems with larger number of nucleons 
with these forces to be incorporated, i.e. ignoring the contribution of AN etc. forces. 
However there is still another possibility to proceed for many-body systems. If the 
contribution of such many-body forces will occur to be significant one can reformulate 
the whole scheme in terms of collective a- (and other) fields interacting with relativistic 
nucleons and dibaryons, i.e. in the spirit of Walecka-Serot model of hadrodynamics. The 
only essential difference in this point from the Serot-Walecka model is adding dibaryon 
components to the nucleonic ones. Moreover, our approach, contrary to the conventional 
OBE-approach, leads inevitably to a strong enhancement of collective cr-field in nuclei 
and thus to validity of Walecka-type model. 

Numerous modern experiments could corroborate these results. In particular, 
according to the recent experiments ^He(e, e'pp) jB] and their theoretical interpretation 
on the basis of fully realistic 3A^ calculations, the cross sections for the ^He(e, e'pp) 
process cannot be explained within a fully realistic 3A^ model incorporating the process 
where a proton knocked out by a virtual 7-ray photon is sequentially rescattered by 
the second nucleon. These calculations show that the leading contribution comes 
from the one-step process where a virtual photon is absorbed by one nucleon in the 
target, while the remaining two nucleons (which are spectators in this process) are 
emitted (in coincidence) at the second step. This important conclusion has been 
further corroborated in recent experiments at the Jefferson Laboratory when the incident 
electron beam energy has been increased up to Ee = 2.2 GeV and 4.4 GeV IH- The data 
of the two different experiments give a clear evidence of very strong short-range A^A^ 
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correlation in the ^He ground state. This correlation still cannot be explained within 
the traditional pattern for the 3A^ system. In addition, our approach has recently been 
partially corroborated [SZ| from the other side by considering a model for 2ti production 
in pp collisions at Ep = 750 and 900 MeV. The authors have found that almost all 
particle energy- and angular correlations (e.g. vr+vr", pp, npp and etc.) can be explained 
quantitatively by assuming that tt+tt" production occurs through the generation of an 
intermediate light cr-meson with the mass m„ ~ 380 MeV and rather narrow width. 
These values generally agree with the parameters adopted in our A^A^ model [011221 and 
drastically disagree with the values assumed in OBE and other potential models. 

Very interesting general implication of the results presented here is their evident 
interrelation to the famous Walecka hadrodynamic model for nuclei jSHj. It is well known 
that the Walecka model describes nuclei and nuclear matter in terms of the scalar a 
and vector tu-fields, where the cx-field gives the attractive contribution, while the vector 
a;-field balances this attraction by short-range repulsion. It is very important that 
both basic fields exist (in the model) as the explicit degrees of freedom (together with 
relativistic nucleons) in contrast to conventional meson-exchange models for nuclear 
forces, where mesons appear as the carriers of forces rather than as the explicit field 
degrees of freedom. Our approach does include the cr-meson (and potentially the u- 
meson) degrees of freedom in an explicit form similarly to the Walecka model. Moreover, 
since the average kinetic energy of the 3A^ system is high in our model (it is higher than 
that in the conventional OBE approach by a factor of more than 2), nucleon motion is 
closer to the relativistic case, and thus the similarity with the Walecka model is even 
closer. 

These general arguments give an additional strong support for the 2N- and 3A^- 
force model presented here. Quite independent numerous arguments in favour of this 
approach in nuclear physics are presented in the subsequent paper 
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